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Usual paradigm in the theory of electron transport is related to the fact that the dielectric per¬ 
mittivity of the insulator is assumed to be constant, no time dispersion. We take into account 
the “slow” polarization dynamics of the dielectric layers in the tunnel barriers in the fluctuating 
electric fields induced by single-electron tunneling events and study transport in the single electron 
transistor (SET). Here “slow” dielectric implies slow compared to the characteristic time scales of 
the SET charging-discharging effects. We show that for strong enough polarizability, such that the 
induced charge on the island is comparable with the elementary charge, the transport properties of 
the SET substantially deviate from the known results of transport theory of SET. In particular, the 
coulomb blockade is more pronounced at finite temperature, the conductance peaks change their 
shape and the current-voltage characteristics show the memory-effect (hysteresis). However, in con¬ 
trast to SETs with ferroelectric tunnel junctions^- 2 here the periodicity of the conductance in the 
gate voltage is not broken, instead the period strongly depends on the polarizability of the gate- 
dielectric. We uncover the fine structure of the hysteresis-effect where the “large” hysteresis loop 
may include a number of “smaller” loops. Also we predict the memory effect in the current-voltage 
characteristics I(V), with I(V) ^ —I(—V). 

PACS numbers: 77.80.-e,72.80.Tm,77.84.Lf 


I. INTRODUCTION 

The single electron transistor (SET) is one of the most 
studied nanosystem,® This is the simplest device where 
strong electron correlations and quantum nature of elec¬ 
tron can be directly observed. It consists of two elec¬ 
trodes known as the drain and the source, connected 
through tunnel junctions to one common electrode with a 
low self-capacitance, known as the island. The electrical 
potential of the island can be tuned by a third electrode, 
known as the gate, capacitively coupled to the island, see 

Fig-0 

For decades there was a paradigm in the theory of elec¬ 
tron transport at the nanoscale related to the fact that 
the dielectric permittivity of nanojunctions was assumed 
to be constant, without any time dispersion. 7 !^ How¬ 
ever, this paradigm is not always true. A number of phys¬ 
ical processes contribute to the polarization of dielectrics. 
Some of them are fast and some are slow compared to the 
time scales of electric field change in the nano junctions. 
Recently, there was a progress in the development of new 
types of dielectric materials with strong and at the same 
time very slow response to the external electric field . 11 ^ 
The SET is a perfect device where this physics can be 
studied. This is related to the fact that the charging¬ 
discharging effects in the SET are controllable and have 
well-defined time scales. 

The Coulomb blockade suppresses the electron trans¬ 
port except for values of the gate voltage where electrons 
sequentially tunnel one by one through SET from source 


to drain. The electric held in the tunnel junctions is 
changing in time while electrons tunnel through the is¬ 
land. The dielectric layers in the tunnel junctions are po¬ 
larized at finite electric held. The usual assumption in the 
theory of SET is related to the fact that the polarization 
of any dielectric layer in the tunnel barrier follows the 
electric held in time: P(t) = a£(t ), where the constant 
a is the dielectric permittivity of the dielectric layer. It 
follows from the last expression that the capacitance C of 
any tunnel junction in the SET is related to the geometric 
capacitance as C = eC^°\ where e = (1+47tq:). And 
this is the only place where the polarization appears in 
the theory of SET. However, these relations have limited 
applicability. In general, the polarization of the dielectric 
is nonlocal in time: P(t) = x(t ~~ T)£(r)dT, where 
x(t) is the dynamical electric permittivity. [Here we as¬ 
sume the linear response regime.] The time dependence 
of function \(t) implies that tuning of dielectric polariza¬ 
tion P(t) by an electric held can not be done arbitrary 
fast. This is happening, for example in dielectric mate¬ 
rials with polarization being due to shift of heavy and 
inert ions. 

The response of polarization P (t) to the external held 
is characterised by the time-scale t p , the decay time of 
function \(t). The second characteristic time-scale in 
the problem: the time of the electric held correlation, t £ . 
For t p -C t £ the polarization has the form P (t) ~ a£(t), 
where a = J^° x(t)c/t. I n the opposite case, r £ <C r P , 
the polarization P(f) does not follow the electric held 
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FIG. 1. (Color online) The equivalent scheme of single elec¬ 
tron transistor (SET). 

S(t) instantaneously and it has the form 

P (*)«<*<£}, (1) 

where {£) is the electric field averaged over the time scale 
t p . It follows from Eq. (|T|) that the simple relation for 
capacitance, C = eC^\ is not valid at shorter times. 
Therefore the theory of single-electron tunneling in the 
SET should be modified and this is the main goal of our 
paper. 

The characteristic time of charge relaxation in the SET 
is t £ = R s C E , where R s is of the order of the bare tunnel 
resistance of the left and right tunnel junctions and C* s 
is the sum of all the capacitances, see Fig. [lj The time 
scale t £ is in the range of dozens of nano- to picoseconds 
depending on the system geometry and materials. The 
switching time of a dielectric material, r P , is in the range 
of seconds to femto-seconds depending on the material 
and the particular physical process behind the polariza¬ 
tion phenomena. 

Therefore the regime of “slow” insulator, t £ -C t>, 
is very important for SET-devices.However, there is 
paradigm that the existing theories with t p -C t £ sat¬ 
isfactory explain most experiments with SETs. What 
is the justification for new theory? The answer is sim¬ 
ple: the effects discussed in this paper are especially pro¬ 
nounced in SETs when on average the polarization of a 
dielectric tunnel junction in the SET is strong enough 
meaning that the charge induced on the grain by the po¬ 
larized dielectric is of the order of the electron charge. 
This condition can be reached for large enough dielectric 
permittivity e only. How large we will discuss below. 

Recently we have found a number of transport effects 
in the SET with slow ferroelectric in the capacitors, see 
Refs. |T]and(21 In particular, we investigated the memory 
effect in this SET. Here we uncover new physical phe¬ 
nomena and report about the memory-effect (hystere¬ 
sis) where conductance periodicity in the gate voltage is 
not broken. Instead, the period strongly depends on the 
polarizability of the gate-dielectric due to the linear de¬ 
pendence of the polarization on the external field in the 
dielectric. Also, we uncover the unusual fine structure of 
the hysteresis-effect, where “large” hysteresis loop may 
include a number of “smaller” loops. We predict that the 
memory effect exists in the current-voltage characteris¬ 
tics, meaning that I{V) ^ —I(—V) for a given memory 
branch even at V g = 0. The last two effects may exist 


in the ferroelectric SET, however non of them have been 
found before. 

The paper is organized as follows. In Sec. [TT] we dis¬ 
cuss the general properties of SET with slow dielectric 
and the methods for investigation of transport proper¬ 
ties. In Sec. m we investigate the SET with slow di¬ 
electric located in the gate electrode at zero bias voltage, 
V2 — V1. In Sec. [TV] we consider the case with slow dielec¬ 
tric in the left and right tunnel barriers of the SET and 
uncover the memory effect in the current-voltage charac¬ 
teristics, I{V). Finally, in Sec. [V] we discuss the validity 
of our approach and the requirements for slow dielectric 
materials which are necessary to observe the effects pre¬ 
dicted in this paper. In the same section we show that 
the Coulomb blockade in SET with slow dielectrics is less 
affected by temperature. 

II. ELECTRON TRANSPORT THROUGH SET 
WITH SLOW TUNNEL BARRIERS 

Here we consider the theory of SET with slow barriers. 
In the following it is convenient to distinguish between 
the geometrical junction capacitances and the low- 
frequency capacitances Ci that include the slow dielectric 
response. The difference between them, aside from the 
unimportant geometrical factor, is 

ACi = Ci- C t (0) = atSi/di, (2) 

where is the dielectric polarizability of the i-th junc¬ 
tion (z = 1 , 2 , < 7 ), Si — the junction surface area and di 
— the effective electrode-island distance. 

We assume that the electrodes are biased with the volt¬ 
ages VT = —V/2, V 2 = V/2 and V g . The grain potential 
(f>(n) at a given number of excess electrons n can be found 
balancing the induced charges: 

ne = £ Cf } (*(«) - V-) + £ AC i «0) - V,), (3) 

i i 

00 

(0) = £ Priori), (4) 

n =—00 

where p n is the probability to find n excess charges on the 
grain. Two terms originate in ([ 3 ]) because we distinguish 
the electric field produced by the capacitance C^ ]> and 
the contribution due to polarized dielectric with slow re¬ 
sponse. So the terms proportional to the coefficient ACj 
in Eq. Q can be considered as charges induced on the 
grain by the polarized dielectric layers that are constant 
in tunneling events. 

The probability distribution p n in the steady state can 
be found using the detailed balance equation^ 

Pn E n ^ n+1 = Pn+lT n+1 ^ n , (5) 

where the rate r n_> ' Il+1 [V, n, (0)] describes the change 
of grain charge from n to n + 1 electrons. The electric 
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current has the form 

OO 

i = e Y, Pn [rr”" 1 - rr n+1 ] • (6) 

71 —— OO 

Here the lower index of T refers to the tunneling rate 
corresponding to the particular tunnel junction, s = 1 or 
2 and the rate T in Eq. ([5| is equal to Id + T 2 . Solving 
Eqs. ([3])-([5]) self-consistently we find the current-voltage 
characteristics of the SET using Eq. <©• 

We use the “orthodox” theory to calculate the 
Coulomb-blockade peaks in the differential conductance 
of the SET. The calculation of T-rates requires the knowl¬ 
edge of the difference in the electrostatic energies when 
the number of excess charges on the grain differ by 
one elementary charge: A= U(n ± 1) — U(n). If 
the polarization in dielectric layers on electron jumps 
follow (f> adiabatically, Pi = aj(</> — Vi)/di, we have 
AU^ = E c ( 1 ± 2 n), where E c = e 2 /2C s with all the 
capacitances C s = Ci being properly renormalized, 

Ci = Cj( 1 + 47rai). However, for slow dielectric layers 
the polarization Pj = aj((tf>) — Vi)/di stays constant dur¬ 
ing the tunneling, and for the energy difference we find 
(see App. 0 

AU±=eW (l±2RT2^P i S i /e) 1 (7) 

where E^ = e 2 /2C^\ d 0) = ^ Cf 5 and PA = 
A Ci({<l>)-Vl). 

The work done by the leads and the gate to transfer 
an electron to/from the grain remains the same as in the 
“orthodox” theory except for the fact that only the geo¬ 
metrical capacitances Cj <r> should be taken into account. 
This implies that for temperature T —> 0 the effective 
ground state free energy is defined as 

F 0 = E ( 0) min(n - Q' /e) 2 , (8) 

n 
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FIG. 2. (Color online) Graphical solution of Eq. (jTo]) show¬ 
ing three possible solutions for an average grain potential (<j)) 
at a given gate voltage V g . Parameters are: Q = — 0.07|e|, 
C ( g 0) = 0.5Cz°\ AC g /Ci 0) = 0.5, and T = 0.4£' o) . The 
three distinct solutions for {(/)) at a given Q o correspond to 
the memory effect instability. 



Q/|e| 


FIG. 3. (Color online) Conductance peaks for A Cg/C^ = 
0, 0.3, 0.6. The “unit” of conductance Gti is the conductance 
of the first tunnel junction of the SET. Parameters are: ca¬ 
pacitances C( 0) = 0.3Ci o) , C^ 0) = 0.5 Ci 0) and C^ 0) = 0.2 
and temperature T = 0.'2Ei°\ The slow dielectric in the gate 
capacitor modifies the shape of the conductance peaks but 
preserves the periodicity in parameter Q in contrast to the 
SET with ferroelectric in the gate capacitor.® 


where the effective gate-induced charge Q' is 

Q'= -C^V g + Y l ^m-V i ). (9) 


Below we use the notation Q = —C g V g for the tradi¬ 
tional gate-induced charge. We show that although the 
effects of slow polarization are far from being a simple 
renormalization of capacitances C^ 1 —> Ci, the conduc¬ 
tance periodicity in Q holds and maintains its period |e| 
for any values of the parameters A Ci. 

The detailed balance equation ([5]) can be solved an¬ 
alytically for the set of voltages V g near the “degener¬ 
acy points”, where the ground state energy of the SET 
changes from n to nil excess charges. The last condition 
requires the effective charge Q' to be close to e(n + 1/2). 
In this case the only two probabilities p n are finite while 
the other probabilities are exponentially suppressed by 
the factor e~ E ° / T . In order to illustrate the origin of 
the memory effect, we will focus on the degeneracy point 
between n = 0 and n = 1 at Vi i2 = 0. Using Eqs. <§-@ 
we find for the average potential (<f>) 


n F [{ 1 - 2 Q'/e)E^] = e(<t>)/2E™ + Q'/e, (10) 


where n F is the Fermi-function. Equation (10) has one 


or three solutions for a given gate voltage Q. The latter 
case is shown in Fig. [2j The presence of three distinct 
solutions for the average potential {</>) at a given param¬ 
eter Q indicates the memory effect instability. Using the 
graphical solution of Eq. ( [To| we estimate the criteria for 
the memory effect instability, ACi/C^ > 2 T/E^\ 
This criterion corresponds to the critical value of AC's 
when the memory effect just appears, see Eq. (32) below 
for the exact expression. 
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AC is A C g . 

For A C g = 0 the conductance is a periodic function of 
the effective gate voltage Q, see the gray curve in Fig. [3] 
The conductance peaks are well fitted by the orthodox 
theory where near the peak maximum the conductance 
is 


a 
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FIG. 4. (Color online) Memory effect instability in the SET 
with slow insulator in the gate capacitor, (a) and (b) the 
conductance branches corresponding to the increasing and de¬ 
creasing parameter Q and A C g /C^ 0> = 0.6, 1, 2, 4, (c) polar¬ 
ization and (d) the average grain potential (arrows show the 
direction of Q evolution for a given branch) for A C g /C^ = 4. 
Grey lines show stable and unstable branches of polariza¬ 
tion and the average potential. Parameters are: capacitances 
C™ = 0.3 C®, C^ 0) = 0.5 C™ and C™ = 0.2C® and tem¬ 
perature T = 0.2_Ec 0) as in Fig. [ 3 ] 


III. SET WITH SLOW INSULATOR IN THE 
GATE CAPACITOR 


A. Numerical study of electron transport through 

SET 


G (0) (dQ (0) ) 


_ eSQW/C^T 

2(R 1 + R 2 ) sinh(e<5O(°)/Gi 0) T)' 


( 11 ) 


Here SQ^/e = min k [-C (0) V g /e - (2 k + l)/2] < 1. 

At finite but small A C g , when the induced charge on 
the island due to polarization is smaller than the elemen¬ 
tary charge, the conductance peaks change their shape, 
but preserve their amplitude and position (see Fig. [3]). 

The opposite case, with dielectric polarization being 
strong enough to induce the charge on the island of 
the order of the elementary charge or larger, is more 
interesting. In this case the conductance peaks show 
the hysteresis and their shape depends on the direction 
of Q-evolution, see Fig. [4] The hysteresis appears for 
A C g > 2T/ E^ (see Eq. 25). Despite the memory 

effect the conductance remains periodic in the renormal¬ 


ized gate voltage Q = —(Cg 0 ' 1 + A C g )V g with the same 
period |e| for any A C g . This behavior is in striking con¬ 
trast to the SET with ferroelectric in the gate where due 
to the nonlinearity of polarization-electric field depen¬ 
dence the periodicity of conductance is broken, see Ref. [21 


Now we discuss the structure of the memory effect. 
Above the critical value of A C g there are many brunch- 
solutions of the self-consistency equation for the average 
grain potential, Eq. Q, for the given temperature, bias 
and gate voltage. The question is - how to choose the 
right branch? Figure [5] provides an answer to this ques¬ 
tion. According to the branching theory ® the j um ps oc¬ 
cur at the “branching points” where the observable has 
an infinite derivative in parameter Q. On the other hand, 
the branch should correspond to the minimum of some 
effective energy functional. In our case (no bias) the role 
of the effective energy plays the free energy 


F = -TlnZ, 


z = exp 

n 


E {0) (n-Q'/e) 2 \ 

T )■ 

( 12 ) 


Here we study electron transport through SET numer¬ 
ically. We consider the SET with slow dielectric layer in 
the gate capacitor. This set-up is the most favourable 
for experiment since in this case there is no electron tun¬ 
neling through the gate electrode and it can be arbitrary 
thick to allow a wide choice of dielectric materials. More¬ 
over, as we will show in the following Sec. |IV[ at V = 0 
by considering the gate capacitor we still preserve all the 
qualitative effects introduced by slow dielectrics in a gen¬ 
eral case. 

Thus, for a time, we assume that the only non-zero 


For zero temperature it reduces to the free energy Fq 
discussed above. 

The plots of the free energy have a similar dependence 
on the parameter Q as the zero-bias conductance G. To 
illustrate this point we show in Fig. [6] the free energy for 
ACg/C^ = 0.6, 1.3. Figure jit shows that the conduc¬ 
tance branch between points “A” and “B” is metastable: 
the free energy for this curve is larger than the free en¬ 
ergy for branch below. However, during the adiabatically 
slow increase of parameter Q the system does not switch 
to the lowest branch at point A, instead it may go up 
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FIG. 5. (Color online) Memory effect: Plots (a)-(d) and (f) show the conductance for A C g /= 1.3, 3.3, 5.3, 10, 20 for stable 
and unstable branches of Eq. 0 for the average grain potential. Plot (e) shows the polarization for A C g /C^ = 10. Arrows 
indicate the position of hysteresis jumps for particular branch with increasing Q. All plots are shown at fixed temperature 
T = 0.2J5c 0) . 


to the metastable branch. The same applies to all other 
plots in Fig. [6j The external perturbation can drive the 
system to outside of the metastable branch before the bi¬ 
furcation point. Usually the role of this “perturbation” 
plays the Langevin forces induced by the thermostat. In 
this case the jumps occur randomly within the same re¬ 
gion before the bifurcation point. This scenario is typical 
for any hysteresis. 

Intuitively one may suppose that if conductance 
“jumps” from one branch to another the final branch 
should have the lowest possible free energy for the pa¬ 
rameter Q corresponding to the jump. Indeed, this is the 
case in Figs. [5ja)-(c) . However, in Figs. E d ) and (f) this 
rule is violated. The system could jump, for example, to 
the point marked by the red-ball in Fig. [5jd) , instead of 
finishing at the point marked by the grey-ball which has a 
larger free energy. However, this energetically favourable 
transition is “forbidden”: while continiously changing the 


polarization in such a process the system would have to 
pass the energy barrier of approximately /4 (free en¬ 
ergy maximum). Thus the higher order jumps (over the 
average charge difference) are suppressed by the factor 
exp(-4 0) /4T). 


B. The fine structure of the memory effect 

Doing numerical studies of memory effect we assumed 
that parameter Q increases (or decreases) monotonically 
from minus to plus infinity (or vice-versa). However, for 
large enough parameter A C g , when polarization induces 
more than one electron on the grain, the hysteresis loop 
depends on the interval where the parameter Q changes. 
This is shown in Fig.[T]with two possible hysteresis loops: 
The red hysteresis loop corresponds to back and forth 
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FIG. 6. (Color online) Free energy in Eq.(T2} for A C g /C^ 0) = 
0.6, 1.3 and temperature T = 0.2The shape of free 
energy plot is similar to the conductance G(Q) plot in Fig.|5k. 
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FIG. 7. (Color online) Memory effect in (a) the conductance 
and (b) the polarization of the gate-insulator. The red hys¬ 
teresis loop corresponds to back and forth change of param¬ 
eter Q in the interval (—2,1), while the blue curve corre¬ 
sponds to (—2,2) interval. Parameters are: A C g = 7.5C^°\ 
T = 0.2 7?c°\ while C^°\ i = 1,2, <7 and Rj, j = 1,2 similar 
to Fig. [3] 


change of parameter Q in the interval (—2,1) while the 
blue curve corresponds to the interval (—2,2). In the 
second case the larger hysteresis loop “includes” smaller 
loops. As a result, the understanding of memory effect 
at finite intervals of parameter Q evolution requires con¬ 
sideration of all branches of the SET observables such as 
conductance and polarization. 


C. Analytical description of the conductance peaks 
and the memory effect 


Here we present the analytical description of transport 
properties of SET. At V = 0 and within the two-state 
approximation the form of the conductance peaks G(Q) 
can be found using Eq. © with the proper substitution 
Q (0) —>■ Q ', where Q' is defined in Eq. ([9]): with this 
substitution we have for conductance G(Q) = G^°\Q'). 
For average potential, generalizing Eq. ([T()|, we obtain 


(<t>) 



-tanh 

2 


^ Ej 0) SQ'^ 



(13) 


where 5Q'/e = min k(Q'/e 
Eqs. (131 with ([9]) we find, 


(k + 1/2)). Combining 


SQ’ 


ci 0) 


e A C g 

2 cF 


tanh 


|^ 4 0) SQ'^j 


SQ, 


(14) 


where SQ = Q — (k + l/2)e is the deviation of parameter 
Q, k is the same as for SQ’ and Gy. = cJj 0 -* + A C g . It 
should be noted that the above equations are valid for 
any SQ as long as SQ' <C 1. 


1. Small polarization 


Here we discuss the limit of small polarization, mean¬ 
ing that the induced charge on the island is small com¬ 
pared to the elementary charge e. Using the small pa¬ 


rameter, A Cg/ 
second order 


-C 1, we expand Eq. (14) up to the 


r (o) 

0 5Q') 0 = 5Q 

Ge 

(<SQ , )i = WOo + l^tanh 


^ 0) (gQ')o ^ 


(15) 

(16) 


The conductance now may be found by substituting 
SQ^ with (JQ^o.i in Eq. [ll] 

G(SQ) = GW(SQ') (17) 

The numerical calculations in Fig. [8])a) show that the 
first order approximation, Eq. |l6| ), well describes the 
peak shape for small parameter A C g /C^ ss 0.1, while 
the zero order approximation is not sufficient. We note 
that parameter A C g /Cg°' > and thus the renormalization 
of the conductance period over V g can be arbitrary in 
this approximation. 


2. Amplitude and form of the conductance peak in the 
hysteresis regime 


Solution of Eq. (14) becomes ambiguous for large val¬ 


ues of parameter A C g , where conductance G(Q) acquires 
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FIG. 8. (Color online) (a) Numerical solution for conductance peak for A C g /= 0.05 (blue line), orange line is the zero 
order solution from Eq. ( |15| ), and the red line is the first order solution from Eq. (161. (b) Numerical and analytical solutions 

~l / f iW) 


of conductance for A C g /= 0.1. The first order approximates well the conductance peak at small A C g . (c) Conductance 
at critical A C g , where hysteresis appears, (d) Conductance hysteresis, (e) Amplitude of conductance peak vs A C g . Points 
represent the numerical solution; the red curve is G ma . x = l/2(J?i + R 2 )', and the orange curve shows Eq. ( |20| ) for Gmax- 
Parameters are: T = 0.06 Ec°\ C.-°\ i = 1,2, <7 and Rj, j = 1,2 similar to Fig. [ 3 ] 


hysteresis. In this case the form of conductance peaks 
becomes nonsymmetric and the conductance G(Q) has 
a maximum at the branching (bifurcation) point corre¬ 
sponding to the jump of the polarization. The bifurcation 
points in Eq. (141 can be found as follows 


d Ln ' CE 

dQ' P Gi 0) 

that reduces to 


cosh 


e A C a 


2d 0) 


tanh 


'4 0) 

T e 


= 0 , (18) 


( 8Q') Y 

T e 


Ei 0) A Cg 
2 T Cy, 


(19) 


The two solutions of Eq. (19) correspond to the in¬ 


creasing and decreasing evolution of parameter Q (solu¬ 
tions with 8Q' < 0 and 8Q' > 0 respectively). These two 
solutions result in mirror-reflected shapes for the peaks, 
so we focus only on the decreasing parameter Q. For 
conductance maximum we find 


arccosh 


f^max — 


E < 0) A Cg 
2 T C s 


2(-Ri + R 2 ) 


I E j 0) A C g ( 

Y 2 T cy V 


E ( c a) ACg 
2 T Cs 


- 1 


( 20 ) 


The predicted conductance maximum amplitude varia¬ 
tion is shown in Fig. [ 8 } One can see that the curve breaks 
at critical value of parameter A C g indicating the start of 
the hysteresis regime. 


We note that since within the scope of the two-state 
approximation and for A C g above the critical value the 
Eq. (20) gives exact maximum, its applicability depends 
only on temperature. At finite A C g the conductance 
maximum does not exactly correspond to a degeneracy 


point 8Q' = 0, but still S(Q') m ax 1 for T <C E‘, 
For example, for temperature T = O.OOFc 0 ^ and AC, 


(0) 


00 we have S(Q') max /|e| ~ 0.1 <C 
consideration is valid (see Fig. [ 8 ]). 


1 , meaning that our 


Now we find the form of conductance peaks. Expand¬ 
ing Eq. (141 up to the second order near 8(Q') ma x we 
obtain 


A 0 + A 2 (8Q'-6(Q') max ) 2 = 8Q, 


( 21 ) 
































where 


where we introduce the parameter 


, eT C E 

0 = Wc^ ar 


( [wac;\ 
\\2T C E J 


and 


eAC g f 2T 

2d 0) V 4 0)A <V 


( 22 ) 


AC's = ^ ACj. (28) 

*=l,2,g 

Here we explicitly show that the functions Q' and (</>) 
depend on voltage V. In general, this dependence results 
in an additional contribution to the conductance propor¬ 
tional to dQ'/dV: 


An — 


Ei 0) C E 
eT 


1 1- 


2 T C s 
EW AC/ 


(23) 


It follows from Eqs. (21) and (17) that the conductance 
derivative in SQ diverges as 1 J\fx near its maximum 
value. 


3. The peak form, at the bifurcation point 

To find the conductance peak at the critical value of 
parameter A C g we expand the hyperbolic tangents in 
Eq. (14) up to the third order. As a result we obtain 


G(Q, V) = 


dlW(Q',V) 
dV ~ 

G^{Q',V) 


dI^(Q',V) dQ’ 


dQ ' 


dV ’ 


(29) 


where ( Q , V) is the current in the orthodox theory, 

generally not limited by the two-state approximation. 
However, the current I is zero for zero bias voltage for 
any Q , therefore the last term can be omitted at V = 0. 
This explains why in two-state approximation we can cal¬ 
culate the conductance by replacing Q by Q' in Eq. (11) 
of the orthodox theory. 


For zero voltage, V = 0, Eq. (27) reduces to 


6Q’ 


AC g ( 4° } 
Cf ] \ 2 T 



^E ( c ° } SQ'^j 

= SQ. (24) 


The linear term equals zero at the critical point. For 
critical polarizability of the gate-insulator we find 

A C { g c) = C^\E^/2T - I) -1 . (25) 

Also we find that 


eT 

SQ' = 


E, 


( 0 ) 


\ 




(26) 


Using Eq. © we find that the peak maximum can be ap¬ 
proximated with the function l/(l + cc 2 / 3 ) (here x ex SQ), 
while the derivative diverges at the conductance maxi¬ 
mum as 1 j yfx. As follows from Fig. ^lc) and Eq. (26) 


this approximation for conductance works well near its 
maximum value only. 


IV. SINGLE ELECTRON TUNNELING 
THROUGH SLOW DIELECTRIC LAYER 


A. Conductance peaks with slow dielectrics in all 
capacitors 

Here we consider the general case, with slow dielectric 
layers in all capacitors with polarizabilities ACi, AC 2 
A C g . Using Eq. (J9| we find 

Q' = Q + AC S x (cf ) }(Q',V)^{AC 2 ~AC 1 )^, (27) 


Q' = Q + AC e x {<f>){Q'). 


(30) 


Then 


SQ' 1 + 


ACe \_ 

Ui 0) J 2 c (°) 


eACs tanh (^)=SQ. 


(31) 

As we can see, the only distinction of the Eq. 131] from 


Eq. 14 is the replacement of AC g with AC's- It follows 
that for V = 0 the SET with slow insulators in tunnel 
junctions behaves qualitatively similar to the only A C g > 
0 that was considered previously. The only difference is 
related to the fact that the slow dielectric in the gate 
capacitor renormalizes the period of the Q-oscillations of 
conductance while slow dielectrics in all other capacitors 
of the SET do not. 

Now we can generalize our results for positive A C g > 0 
obtained earlier. In particular, the critical polarization, 
where memory effect in the conductance G(Q) first ap¬ 
pears, becomes the integral quantity, see Eq. (28), that 
includes properties of all the slow dielectric layers: 


AC^ ] = C^\Ef>/2T- 1) 


-1 


(32) 


The amplitude of conductance peaks can be found using 


the substitution, A C„ —5> ACe in Eq. (20). The shape of 


the peaks can be obta ined us ing the same substitution in 
the equations of Sec. IIIC3 where still SQ = —(C g °' > + 
A Cg)Vg. 


B. Memory effect in current-voltage characteristics 

Above we discussed the properties of SET with slow 
dielectric barriers, related to the variation of the gate 
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FIG. 9. (Color online) Current-volt age characteristics, I(V) of SET with zero gate-capacitance, C g = 0. Plot (a) shows I(V) for 
ACi/C^ 01 = 0.0, 0.25, 0.5, 0.75,1.0,1.25. Smoother curves correspond to smaller AC/. Plot (b) corresponds to AC// C^ = 1.5. 
The jumps in /(V) in (b),(c), and (e),(f) correspond to the memory effect: the branch depends on the direction of voltage 
change. Plot (c) shows I(V) for AC/ / Cy 0> = 3.5. The /(F)-curve can have many hysteresis loops depending on the amount 
of electron charge induced on the grain by the dielectric polarization. Inserts in (b),(c) show the details of the hysteresis. Plot 
(d) shows I{V) for ACi/C^ = 1.25, A Ci/C^ = 0 (black curve) and AC//C .^ = ACi/C'^ 0 ‘ > = 1.25 (orange curve), while 
plot (e) shows the graphs for ACi/C® = 1.5, AC 2 /= 0 (black curve) and ACi/C^ = AC 2 /= 1.5 (orange and blue 
curves). Plot (f) shows I(V) for ACi /= ACb/C^ 0 ' = 3.5. Parameters are: T = 0.06 i?c 0 \ C{ 0) = 0.6Ci; 0 \ C^ 0 ' = 0.4Cs°\ 
and Ri, i = 1,2 similar to Fig. [ 3 ] The unit of voltage is i?c°V|e|, and the current is normalised to /\c\Rti. 


voltage V g at bias V = 0. In this subsection we instead 
concentrate on the current-voltage characteristic I(V) of 
SET in the case of electron tunneling through slow insu¬ 
lator in the left and the right capacitors, see Fig. |T) We 
neglect the gate to simplify the situation, putting thus 
C g = 0. Such systems have been extensively studied in 
experiments over the last two decades. They can exhibit 
coulomb blockade at room tenxperature dTTl ITTHIS]) and 
their ease of fabrication makes a wide range of barrier 
materials available for experiments. Following Ref. [TT; 
we consider the current-voltage characteristic of the SET 
in a wide range of bias voltages. 

The typical current-voltage characteristics I(V) are 
shown in Fig. [9] in Fig. |9])a)-(c) the coefficients AC 2 = 0 
and AC/ are finite. It follows that there is a memory 
effect in I(V) at large enough AC/ and this effect de¬ 
pends on the direction of the bias voltage evolution. The 


jumps in plots (b) correspond to the regions of hysteresis 
while the arrows show the evolution of voltage. Plot (c) 
shows the hysteresis in I(V) for ACi/C^ = 3.5. The 
current-voltage characteristics may have many hystere¬ 
sis loops, depending on the amount of electron charge 
that the dielectric polarization may induce on the grain. 
The hysteresis in the current-voltage characteristics ap¬ 
pears for the first time for parameter AC/ being larger 
than C l i°' > . This is the first critical value of polarization. 
For AC 1 > 2Ci°' ) the second hysteresis loop appears in 
I(V). Therefore this is the second critical value of AC/. 
For larger values of AC/ we expect further increase in 
the number of hysteresis loops. 

Two cases of current-voltage characteristics are com¬ 
pared in plots (e)-(d) : i) finite AC/ and zero AC -2 and 
ii) AC/ = AC/. In both cases the set of critical values 
of AC is the same and for large bias voltage the current- 
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voltage characteristics asymptotically coincide. 

Figure [9] shows that the current-voltage characteristics 
of the SET strongly depend on the direction of bias volt¬ 
age F. Moreover, for a given hysteresis branch 

I(V) ± -I(-V) (33) 

that happens in the absence of Q, notably different from 
the result for a regular SET. 


C. Influence on coulomb ladder 


By coulomb ladder in this section we mean a step-like 
behavior of /(F) in the regime of coulomb blockade. The 
coulomb ladder is often used as an indication of coulomb 
blockade (Ref. [TTJ [TU 177111111) . In the following we show 
how the slow polarization influences the shape of the lad¬ 
der. Again we take C g = 0 and, consider the conditions 
when the ladder is the most pronounced, i.e. T = 0 and 
strongly asymmetric barriers R\ R 2 . At zero temper¬ 
ature tunneling may occur only in the direction of chem¬ 
ical potential drop, that is from the 1-st electrode to the 
2-nd assuming V > 0. Due to the relatively high tunnel¬ 
ing rate through the 2-nd electrode, the number of excess 
electrons on the island is almost always stays at the min¬ 
imum energetically allowed number n m i n . n m i n can be 
determined as the lowest n for which AF^ +1_> " < 0 is 
true, since Af’ 1 n_> " +1 < 0 holds for any n < 0. For a 
given n m i n the current can be calculated as 


I = -A F, n 

eR\ 


»+i 


(34) 


where AR is the free energy change on tunneling through 
the 1-st electrode. For a conventional SET the above 
formula leads to a ladder-shaped /(F) characteristic with 
the step width 


AF s tep = e /cj 0) , 

(35) 

jumps of the current between the steps 


A/ st e P = |e|/Rrd 05 , 

(36) 

and the /(F) slope between the jumps 


dl/dV = C^ 0) F/ci 0) /?i 

(37) 


Introducing slow dielectric into the tunnel junctions re¬ 
sult in some new effects (for the details of calculations 
see Appendix]^]). At V > \e\/C^ slow polarization leads 
to the rescaling of the ladder that may be described by 
substituting the capacitances in Eqs. 35]|37| with the new 
values Ci = C j 0 - 1 + AC). exactly as when dealing with 
a conventional fast dielectric (see Fig. |l0](b)). But con¬ 
trary to the fast dielectric, the slow one shifts the lad¬ 
der, making it asymmetric and, moreover, dependent on 
the direction of the evolution of F, as illustrated at the 

Fig. pM). 


Interestingly, the shift of the /(F) curve in experiments 
is a well-known effect. It is usually accounted for by 
assuming the presence of some additional spurious charge 
Q, induced on the grain (as in Ref. fill and [75)) . However 
the shift that we predict is notably different at least in 
one aspect — it reverses its sign with the direction of the 
evolution of V. 

We stress that the described rescaling and shift of /(F) 
takes place only under specific conditions V > |e|/c|°^ 
and R\ i? 2 - If R 2 are of the same order the intro¬ 
duction of slow dielectric may change the ladder steps 
in a more complex way. Such a situation is shown in 
Fig. [l0](b) where the ladder period do not correspond to 
the one we would expect from the simple capacitance- 
renormalization consideration. If Ri/R 2 is even closer 
to unity, the slow dielectric barriers qualitatively change 
the current-volt age curve as was discussed in the previous 
section (see Fig. [9]). 


V. DISCUSSION 

A. Requirements for dielectric materials 

Here we discuss several possible dielectric materials 
which can be considered as slow insulators. At finite 
external electric field the localized electric charges are 
shifted and the dielectric material is polarized. There 
are several physical processes contributing to the polar¬ 
ization: 1) the shift and deformation of electron-cloud, 
2) the shift of ions in the lattice, and 3) the molecular 
and/or macro dipole reorientation. Electrons, ions, and 
dipoles can form a different polarization. The slowest po¬ 
larization formation corresponds to the electrocalorical 
and migration (electron, ion or dipole) mechanisms with 
the characteristic dispersion frequency being in the range 
10“ 4 — 10“ 1 Hz and 10“ 3 — 10 3 Hz, respectively at tem¬ 
perature T = 300/F The electromechanical mechanism 
corresponds to frequencies 10 5 — 10 8 Hz, while thermal 
mechanism correspond to 10 5 — 10 10 Hz. The dielectrics 
where thermal mechanism is the largest are promising for 
applications in nanostructures and can be considered as 
“slow” dielectrics. 

Dithiol self-assemble monolayers (SAMs) have a static 
dielectric permittivity e(u> = 0) ~ 3 and the character¬ 
istic relaxation frequency ~ 10 4 Hzf 12 These materials 
are good candidates for slow dielectrics. Such dielectric 
layers have been used in double junction SET^!]. The 
hysteresis have not been observed in these experiments, 
but there was a considerable discrepancy between the the 
values of capacitances obtained from the fit of the exper¬ 
imental data with the orthodox model and the ad-initio 
calculations. 

Another promising materials to observe the hysteresis 
are polar crystal dielectrics e.g., BaTiC >3 or KDP with 
static dielectric permittivity e(w = 0) ~ 10 3 and the 
typical relaxation frequency w c ~ 10 6 Hz. 
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FIG. 10. (Color online) Current-voltage characteristics, I{V), of SET demonstrating the coulomb ladder at T = 0, C g = 0. (a) 
I(V) in the regime when the scaling of the coulomb ladder steps at large V is the same for slow and fast dielectric response. 
Here R 2 /R 1 = 10 —3 , C ^ = 0.9 c£\ C^ = 0.1 C^°\ (b) The regime when the periods of coulomb ladders for slow (blue curve) 
and fast (red curve) dielectric of the same static polarizabilities are different. R 2 /R 1 = 0.25, = 0.9 C^, Ci°^ = 0.1 C^\ 

(c, d) The shifts of the ladders arising from the hysteretic behavior of I(V) for SET with slow dielectrics. Arrows indicate 
the directions of voltage change for each curve. Here R 2 /R 1 = 10 —3 , Cj; 0 = 0 .7C^\ = 0.3C^K The unit of voltage is 

\e\/Cg\ current is measured in \e\/C^\R i + R 2 ) 


B. Fast capacitances D. Temperature dependence of the 

coulomb-blockade effects 


Here we discuss the geometric capacitance Cj°\ i = 
1,2 ,g. We assumed that these capacitances has an 
electrostatic origin. However, in rigorous analysis 
they include the high frequency dielectric permittivity 
Coo (usually between 1 and 10). Thus in our consideration 
the slow polarizability on is the difference between the low 
and the high frequency on- As an example, for BaTiOa 
the difference between the high and low frequency per¬ 
mittivities e is ~ 10 3 . This difference is large enough. 


A well-known consequence from the orthodox theory 
of SET is that in order to experimentally observe the 
coulomb-blockade phenomena, the temperature of the 
system should be lower than E c = e 2 /2 C^. Here the 
total capacitance Cs includes dielectric susceptibility of 
the barrier media. In contrast, our numerical calculations 
show that if the dielectric response is sufficiently slow, 
only the ratio E^ /T should be taken into account when 
considering the blurring of the coulomb effects due to fi¬ 
nite temperature. This must result in more pronounced 
blockade for a system with slow dielectric at a given tem¬ 
perature and electrode geometry, as illustrated in Fig. El 


C. Critical polarization 

VI. CONCLUSIONS 


The effects of slow polarization are governed by the 
ratio of ’’slow” and ’’fast” capacitances A Cs/ci°\ If a 
capacitor is fully filled with a dielectric with permittivity 
e(w) than AC^/C^ = (e(0) — e(oo))/e(oo). It follows 

i(0) 


from Sec. 


Ill 


and 


IV 


that at AC's / ~ 1 the strong 

influence of slow polarization may be observed, thus re¬ 
quiring e(0) > 2e(oo). 

The latter requirement become even less strict at 
lower temperatures. In particular, the critical value of 
e(0) (c) A( 00 ) to observe the breakdown of conductance 
peaks goes to 1 as T —► 0 (see Eq. (32)). For the condi¬ 
tions as at the Fig. [she) e(0)^ c ^ « 1.14 e(oo) 


We showed that the dielectric materials at the 
nanoscale demonstrate new physical phenomena. As an 
example we studied the single-electron transistor. We 
found the memory effect in the conductance-gate voltage 
dependence and in the current-voltage characteristics of 
the SET. We uncovered the complex fine structure of the 
hysteresis-effect, where the “large” hysteresis loop may 
include a number of “smaller” loops. We also found, that 
in order to estimate the influence of temperature on the 
electronic transport one should compare T with e 2 /2 
where in C^ the slow part of the dielectric function is 
not included. 
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FIG. 11. (Color online). Temperature dependence of the 
coulomb ladder in I(V) characteristics of (a) regular SET and 
(b) SET with slow dielectrics in the tunnel barrier. In plot 
(b) the SET parameters are c)°' ) = 0.7 C^\ C^ = 0.3 C^°\ 
ACi = 1.4 d°\ AC 2 = 0. In plot (a) Ci = 2.1C^ ) and 
C 2 = 0.3 C^. For both plots R 2 /R 1 = 10 -3 and in (b) 
I(V) is shown for increasing voltage V. The unit of voltage 
is |e| /C^°\ the current is measured in |e|/Cs°^(^i + R 2 ) 
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Appendix A: Calculation of the coulomb energy 
change on electron jumps 


Here we show, how the energy changes A[/„ are cal¬ 
culated. If the number of electrons on the island changes 
from n to n ± 1 in some process, than the electrostatic 


energy change is 


Art 



^(</>- Vi)dqi = 





(Al) 


where qi are the charges of the capacitors and Pi are 
dielectric polarizations in barriers. For fast and slow di¬ 
electrics Pi behave differently during the process of elec¬ 
tron jump. If dielectric response is fast Pi follows (j) that 
results in capacitance renormalization. For slow dielec¬ 
tric layers the polarization cannot change on the electron 
jump timescale and thus dP = 0 yielding 


Art = - 


4>(n± 1 ) 


4>{n) 


(A2) 


4>{n) are calculated using the charge balance equation © 


<t>{n) = 


a 


( 0 ) 


e n 




■ (A3) 


Here Pi are constant and do not depend on n. By insert¬ 
ing Eq. (A3) into (A2) we obtain Eq. Q. 


Appendix B: The shape of the coulomb ladder 


At zero temperature the tunneling rates for the elec¬ 
tron to and from the island are 


1,2 


*n± 1 


e 2 i?i, 2 


(-A F^ n±1 ) 


0 (—AF"^ >n±1 ) , 


(Bl) 


where n is the number of excess electrons on the island 
and tunneling happens through the 1-st or the 2-nd elec¬ 
trode. Free energy changes AFj, 2 on jumps are 


A F™~ 

■4-nil _ e 

/o(0) 

°E 

(i ±(ne 

-Q')±cf ) V^j 

(B2) 

a Fr 

►n±l _ e | 

/o(0) 

°E 

(i ±(ne - 

-q') t c[ 0) v) . 

(B3) 


Consider V > 0. It follows from |B1| that tunneling 
occurs if for some n simultaneously A n+1 < 0 and 
A F™+ l ^ n < 0 (there is no backward tunneling at T = 0). 
These conditions may be combined into 


Q'/e- 1/2 + Cf ] V/e <n< Q'/e—1/2 — V/e. (B4) 


Since the tunneling from the 1-st electrode to the island 
is much slower than from the island to the 2-nd electrode 
(i?i R 2 ), the number of electrons on the island almost 
constantly stays at it’s lowest energetically allowed value 
rimin- The current is, than, 


-pn—m+l-pn+1—m 
1 1 _ ^2 _ 

pn->n+l . pn+l-ln 
1 1 ' 1 2 




(B5) 


1= -e 
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The rest is to calculate n m i n . Since we neglect C g only 
the charge induced by the slow polarization gives rise to 

O' 

_ AC S _ „ , ACrCf } - A C 2 C[ 0) T , 

Q — ^min^ “I" ^ • (B6) 

Os Os 

n-min can be determined from the equation 


1 


2 1 + ACs/ Ce (1 + ACs/Q 03 ) 



where [a:] denote the lowest integer bigger than x. It 


worth noting that the equation B7 predicts multiple so¬ 
lutions for n m i n at V close to the current jump points if 
ACs > 0 (see Fig. [Toj(d)). 

The calculation of I yields 


I(V) = 


1 


e Cs 


fliCs \ 2 d 0 ) 


+ ^mine + C2V . 


(B8) 


The latter formula demonstrates the full renormalization 
of capacitances and a shift in the /( V) as is illustrated 
at the Fig. [lOja). 
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